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We prove that static, spherically symmetric, asymptotically flat soli- 
ton and black hole solutions of the Einstein- Yang-Mills equations are 
unstable for arbitrary gauge groups, at least for the "generic" case. 
This conclusion is derived without explicit knowledge of the possible 
equilibrium solutions. 
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1 Introduction 

In several recent papers H, 0, pi E| , we have studied important aspects of the 
Einstein- Yang-Mills (EYM) system for arbitrary gauge groups. In particular, 
we investigated the classification and properties of spherically symmetric EYM 
solitons (magnetic structure, Chern-Simons numbers) and a generalization of 
the Birkhoff theorem for the non-Abelian case. We also worked out the general- 
ization of the first law of black hole physics (Bardeen-Carter-Hawking formula), 
allowing for additional Higgs and dilaton fields g, || . For other studies of these 
and related topics we refer to J7], g, g, [hJ. 

In the present paper, we prove that static, spherically symmetric, asymptot- 
ically flat solutions of the EYM equations are unstable for any gauge group, if 
they are "generic" (defined in Sec. 2). In a recent letter pj , we have already 
sketched how we arrived at this result for solitons. Here we present details of the 
proof and extend it to black holes. We discuss also some further mathematical 
issues involved. 

This general instability was expected since the Bartnik-McKinnon solutions 
L2| for the gauge group SU(2) and the related black hole solutions 13 [14], 



[15| are unstable |IJ, |l7], [T|, ^9[. A mathematical proof of this expectation 
presents, however, quite a challenge, since one cannot rely on any knowledge of 
the possible solutions (apart from regularity and boundary conditions). 

Our strategy is based on the study of the pulsation equations, describing 
linear radial perturbations of the equilibrium solutions and involves the follow- 
ing main steps: First we show, that the frequency spectrum of a class of radial 
perturbations is determined by a coupled system of radial, respectively one- 
dimensional "Schrodinger equations". Negative parts in the spectrum of the 
effective Hamiltonian imply linear instability. With the help of suitably con- 
structed trial functions, it is then proven, that the spectrum contains always a 
negative part (for "generic" solutions). 

We have recently used a similar procedure to establish the instability of the 
gravitating, regular sphaleron solutions of the SU(2) EYM-Higgs system with 
a SU(2) Higgs doublet (g0|, which have been constructed numerically in fl2lf . 
Our results contain, as a special case, the conclusion of Ref. [£2| for the gauge 
group SU(2). Here, we analyze the regular SU(2) case further. We show that 
the effective Hamiltonian for "sphaleron-like" perturbations has the form of a 
"deuteron" Hamiltonian. 

The paper is organized as follows: In Sec. 2 we recall some basic facts and 
equations of our previous work [0, B , which will be needed in the present anal- 
ysis. In Sec. 3 we then derive the linearized perturbation equations for solitons 
and black holes and bring them into a convenient, partially decoupled form. 
The resulting eigenvalue problem is discussed in Sec. 4 and in Sec. 5 we show 
the existence of unstable perturbations. The "deuteron" interpretation for the 
unstable modes of a SU(2) soliton is presented in Sec. 6. In the appendix, we 
elaborate on mathematical issues, related to the self-adjointness of the effective 



Instability of EYM Solutions 



Hamiltonian and the connection between the negative part in its spectrum and 
unstable solutions of the perturbation equations. 

2 Spherically symmetric EYM fields 

We begin with a convenient description of gauge fields with spherical symmetry 
(for derivations see ||]). 

Let us fix a maximal torus T of the gauge group G with corresponding 
integral lattice / ( = kernel of the exponential map restricted to the Lie algebra 
LT of the torus T). In addition, we choose a basis S of the root system R of 
real roots. The corresponding fundamental Weyl chamber 

K{S) = {H eLT\ a{H) > for all a e S} (2.1) 

plays an important role in what follows. 

To a spherically symmetric gauge field there belongs a canonical element 
H\ € InK(S), which characterizes the corresponding principal bundle P(M, G) 
over the spacetime manifold M, admitting a SU(2) action. If the configuration 
is also regular at the origin, H\ is restricted to a small, finite subset of IC\ K (S), 
which is described in 1. In the present discussion, we exclude (for technical 
reasons) the possibility that H\ lies on the boundary of the fundamental Weyl 
chamber. The term "generic" always refers to fields, for which the classifying 
element H\ is contained in the open Weyl chamber K(S). 

The SU(2) action on P(M, G) by bundle automorphisms induces an action 
on the base manifold M. A SU(2) invariant connection in P(M, G) defines 
an invariant connection in each subbundle over a single orbit of the action on 
M. By Wangs theorem, the induced connections are described by a linear map 
A: LSU(2) — » LG, which depends locally smoothly on the orbit and satisfies 

A 1 = [A 2 ,A 3 ], A 2 = [A 3 ,A 1 ], A 3 = -H x /4n, (2.2) 

where A& := A(r^) and 2iT)~ are the Pauli matrices. These equations imply that 
A + := Ai + iA2 lies in the following direct sum of root spaces L a of LGq'- 

A+ € {^0 L a , 

aes(A) (2.3) 

S(\):={aeR+\a(H x ) = 2}, 

where R+ denotes the set of positive roots in R (relative to the basis S). In the 
generic case S(X) turns out to be a basis of a root system contained in R (see 
appendix A of Ref. j|). 

The LG- valued functions A± on the orbit space determine part of the con- 
nection on P(M, G). Before we give a parametrization of the YM potential in 
a convenient gauge, we fix our conventions in parametrizing the Lorentz metric 
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on M and introduce some further notation. We use standard Schwarzschild-like 
coordinates and set 

ds 2 = -NS 2 dt 2 + A" 1 *- 2 + r 2 (dd 2 + sin 2 tf dip 2 ) , (2.4) 

where the metric functions N = : 1 — 2m/ r and S depend only on r and t. 

A suitably normalized Ad(G)-invariant scalar product on LG will be denoted 
by (•,•)• We use the same symbol for the hermetian extension to LGc (linear 
in the second argument), and | • | means the corresponding norm. Note that the 
original Ad(G)-invariance extends on LGc to 

(X,[Z,Y]) + ([c(Z),X],Y) = 0, (2.5) 

where c is the conjugation in LGc- 

In Ref . S , it is shown the gauge potential A can be chosen to have the form 

A = A + A, (2.6) 

with 

A = k 2 d$ + (A 3 costf- Kxsm.d)d(p (2.7) 

and 

A = Adt + Bdr, (2.8) 

where A and B commute with H\ (i.e. with A 3 ). Since H\ is assumed to 
be generic, its centralizer is the infinitesimal torus LT. Hence, A and B are 
LT-valued and A is thus abelian. 

For the example of the gauge group SU(2), H\ is an integer multiple of 



47r T3: H\ = Ank T3 with k e Z, and the only solutions of (2.2) are Ai = A2 = 0, 
A3 = k T3 and 

Al=WTl+WT 2 , A 2 = Tw T\ ± w r 2 , A 3 = ±t 3 . (2.9) 

The gauge potential A contains a "trivial", abelian part, which decouples 
from the EYM equations. To demonstrate this, let us first construct a convenient 
decomposition of LT. For a given potential we restrict the sum in Eq. (|7|) to 
the smallest subset S of 5(A) for which 

A+e L a . (2.10) 

QGSCS(A) 

Since every rootspace L a is Ad(T)-invariant and since the residual gauge group 
of the potential A is just the torus T, the subset S is unique and depends only 
on the invariant connection. With the help of S we now split LT: 

LT=(£)e(£)\ (2.11) 
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where (S) denotes the linear span of £. The decomposition ( p. 11 ) is indepen- 
dent of the chosen Ad(G)-invariant scalar product B and satisfies 



[<£)\A + ]=0 



This property motivates to set 



A = a + A, 

B = b + B, 



(2.12) 



(2.13) 



A, 



A, 



A 



3 II 



with a, b, Asx € ( £ ) ± and A, B, A 3 || e ( £ ). For our instability proof we adopt 
the following (mixed) gauge: 



.4 = , 6=0, 



lim a = 



(2.14) 



If we now insert the parametrizations (2.4), (2.6) - (2.£), ( 2.13| ), (2.14) into 
the EYM equations, we obtain a system of partial differential equations for the 
metric functions N, S and the YM amplitudes A±, B. As noted above and as 
Eq. (2.12) indicates, the equation for a decouples. Specializing the results of H 
(and using slightly different notation), they read as follows: 

The Einstein equations give two constraint equations for the r derivative 
(denoted by a dash) and the t derivative (denoted by a dot) of m 



i' = ^{NG + r 2 Pe }, 



m = -NH , 
2 



(ft :— 8nG), and the (rr)-equation reduces to 

S' K 



s 



-G, 



where 



G = \\ (NS)- 2 \\ + \ 2 -\A\ +[RA H 



H = Rc(A+, A' + + [B,A+}) 
1 



po 



2r 4 



l^ii 



w 



p± 



\Qi 



(2.15) 

(2.16) 

(2.17) 
(2.18) 
(2.19) 



with 



and 



A = -[A + ,A_]-A 3 



P± = A 3± 



^? B 



"T a 



(2.20) 

(2.21) 
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The YM equations decompose into 



2 fr 2 



NS\S 



B 



-(- 

S\NS 



A, 



A+ , A'_ + [B, A_] + A_ , A' + + [B, A+] = , (2.22) 



±(iVs{AV + [B,A + ]} 



iV 



£,A' + + [£,A H 



[^||,A + ] 



^l-HI 



21 ?* 



NS\ 



[A + ,A_] + [A_,A+] 



0, (2.23) 
(2.24) 



The abelian electric part of the potential satisfies 

^.2 



constant (g (£) ) 



(2.25) 



and hence decouples. 

Eqn. ( 2.24 ) is the Gauss constraint. For static solutions all time derivatives 
disappear, B can be gauged away and the basic equations simplify considerably. 
(For the Bartnik-McKinnon solutions A is of the form (2.S) with w — 0, A311 = 
A 3 = r 3 and A = in (O).) 



3 Perturbation equations 

In this section we study time-dependent perturbations of a given static, asymp- 
totically flat solution of the coupled EYM equations ( E.15|), (J2.16Q , d2.22j ) - 
( |2.25| ) . Regular solutions are "purely magnetic" (A — in (|2.6| )) with vanish- 
ing YM charge (P±=Q±=0 and limj-^oo j°^=0). Unfortunately, this is not yet 
proven with satisfactory weak fall-off conditions, but there is strong evidence 
for this (see jl], g3| for partial results.) The perturbation equations we derive 
hold also for black holes, if their gauge potentials A have the form 



A 



with 



a(r) = Qj_ 



adt + A 

S 



ir 



dy 



(3.1) 



(3.2) 



for a constant vector Q± in (£) (i.e. A = B = in Eq. (2.13)). We call such 
gauge fields "essentially magnetic" . 

From now on A±, N, S, etc. refer to an essentially magnetic equilibrium 
solution and time-dependent perturbations are denoted by SA±,SB, etc.. All 
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basic equations are linearized around the equilibrium solution. In order to 
decouple the perturbation 5a, we impose the additional constraint 5Q±=Q. 



First, we linearize the right hand sides of the Einstein equations (2.15) and 
( |2.16| ). Since B and A± vanish for the equilibrium solution, the first order 
variation of the source G is 



5G = Re ( A' + , SA' + ) - Re ( A' + , [A+,5B]] 



(3.3) 



Here, the last term vanishes, because the property ( pL5| ) of the scalar product 
implies 

-2Re{A' + ,[A + ,5B]) = {[A + ,A'_] + [A_,A' + ],5B), (3.4) 



and the YM equation ( 2.22 ) for the equilibrium solution shows that 

[A + ,A'_] + [A_,AV]=0. (3.5) 

Thus, 



SG = Re ( A' + , SA' + ) . 
The only first order variation for pg comes from 5\!F\\ \ - 

5f ll = ^{A + ,dA_}-~{A_,5A + ] 



(3.6) 



2(^,1, 8F\\). Usine 



"II . "^I 



(see Eq. ( [2.20 )), we have 



Sp e = —Re(i[T h A + ], <5A H 



(3.7) 



(3.8) 



Now we can work out the variation of the first Einstein equation in (2.15). 
With (|3.6|), (3.8) and (2.16) for the equilibrium solution, we find 



5m' = -^5m+~lNRe(A' +> 5A' + )+ Rc( -l[i||,A+], SA+)\ . (3.9) 

For the commutator in the last term we use the unperturbed YM equation 
(HI), i-e. 

4[^||,A+]=A^a;+^'A; + A^, (3.10) 



S' 



whence 

S' 
5m' = — — 5m 



or 



+ f { f ^Re( A' + , 5A + ) \ + | ^NRe(A' + , 8A + ) | (3.11) 
(5mSy = S.^NSRe(A' + ,5A + )\ . (3.12) 
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Therefore, 8m must be of the form 



8m = ^iVRc(A' + , <5A+) + ^r > (3.13) 

where f(t) is a function of t alone. This function is determined by considering 



the variation of the second Einstein equation in (2.15), which reads 



Thus, we have also 



8m = ^NRe{A' + , SA + ) . (3.14) 



Sm=^NRe(A' + ,SA + )+g{r) , (3.15) 



8(^)=^NRe(A' +> 8A' + ). (3.17) 



with a function g(r) of r alone. By comparing (3.13) and ( 3.15| ), we arrive at 
the remarkably simple result 

Sm = ~NRe(A' + , 5A+) , (3.16) 

which generalizes an observation already made in |lq ]. 

The variation of the Einstein equation ( 2.16| ) is immediately obtained with 

(0) 

'S'\ _ K 

Before also linearizing the YM equations, we introduce a suitable decompo- 
sition of A + and <5A + . To do so, we choose a base element e a of the root spaces 
L a and expand the unperturbed A + as well as its perturbation 8A + : 

A+ = ^u; Q e Q , SA+ = ^2 8w a e a . (3.18) 

Then, we have 

6A± = SX± ± i6Y± (3.19) 

with _ ___ 

8x+ = j2 Re ( 5w<x ) e « ' 6Y + = Yl Im ( 5wa ) ea ( 3 - 2 °) 

and the corresponding expansion for 8X_ and 8Y_ with e a replaced by c(e a ) G 
L_ Q , because £A_ = c{8A + ) and thus 

8X^ = c(<5X+ ) , 8Y_ = c{SY + ) . (3.21) 

We call 8X± , 8Y± the "real" (or "gravitational" ) and "imaginary" (or "sphalc- 
ron-like" ) parts of the perturbations SA± . It was shown in B that the unper- 
turbed A + can be chosen to have only a real part. 
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This decomposition will lead to a significant decoupling of the perturbation 



equations. Note in particular, that the variations 5m and Spe in (3.8) and (3.16) 
depend only on the real part 5X + : 



5m = -N(A' + ,5X+), 



Spe 



i[F h A+],6Xi 



(3.22) 
(3.23) 



We consider now the first variation of the YM equation (2.23). Its decom- 
position into real and imaginary parts yields 



' 6X+ = -N5X'l- WLSX' + - 4[A+,^||] 



NS 2 



r[^||,<JA H 



_ . ,VA". - ,i['^ZY\'_ 



and 



NS 2 



■5Y, 



-N 



(Nsy 
s 



s 



[A+,5B]' + i[A'5B] 



(3.24) 



5Y[ + i [A+JB] + - [F h 8Y+] .(3.25) 



The third term on the right hand side of (3.24) is indeed real and can be written, 
using (fO|), as 

\ [A+, 5T\\ ] = \ ad(A + ) ad(A_) 5X+ . (3.26) 



Equation (3.24) can be simplified further. From ( |3.22| ) and the equilibrium 
equation (3.10), we deduce 



- 5NA" 



5m A'{ 



= kN Re ( A' + , 5X+ ) A 



= k Re ( A+ , 5X^ 



(NS)' 
S~ 



A, 



and the Einstein equations (2.15), (p.K:) give 



(Nsy 
s 



r5m + fvr < 



r2^lh A +] 



(3.27) 
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If we use here ( 3.2 2| ) and (3.23), we see that the last two terms in (3.24) can be 
expressed as follows: 






NS 2 



-AW{<^ + ^} W , 



sx 4 



2 2 1 

A' + K fL ( [iJP,, , A+] , SX + ) + [#|| , A + ] K !y < A' + ,6X+) , (3.28) 



where 



NS 2 



M 



(3.29) 



Inserting these expressions into (3.24) gives the following pulsation equation for 
the real amplitude 5X + of the YM field: 



SX+ + U XX 5X+ = Q, 
where the operator Uxx is given by 

Uxx = p 2 + A* 2 ad(^||) - /i 2 ad(A+)ad(A_) 

t' +K /i 2 {l- K rV}(AV, ■ ) 



(3.30) 



A 



2 2 

+ A^ K ^-( [#| h A+], • )+ [i£j|,A+]K^-« 



A 



+ ■ 



(3.31) 



and p* denotes the differential operator 



-iJV5 



9? 



(3.32) 



It is remarkable that the perturbations 5Y± and SB do not appear in (3.30) 
and that the back reaction of gravitation on 8X + can be described by an effective 



potential (last three terms in (3.31)). 

Equation ( |3.25| ) can easily be brought into the form 

5Y+ + U YY 5Y+ + Uyb^/Nt SB = . 



where 



Uyy — p 2 + ^ 2 ad(«.F||), 

Uyb = p*^ad(A + ) + /xad(p*A + ) 



(3.33) 

(3.34) 
(3.35) 



We have thus achieved a partial decoupling, because neither SX + , nor the 
metric perturbations, appear in (3.33J). 
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We proceed with the linearization of the YM equation (2.22). The variation 
of the last two terms is 

A+, [A-,SB] + [A+,SA'_] - [A'_,<5A + ] + conjugate, (3.36) 

which leads (with SA± = 5X± ± iSY± ) to 

A+ , [A_, SB] 1 + i [A+, 6YL] + i [A'_,6Y + 



+ 1 [A+,SX'_] - [A'_,5X+] 1 + conjugate. 



Here, the terms in the first curly bracket are in LT, while those in the second 
are in iLT. The latter are compensated by their conjugates and we find 



with 



Nr SB + UbbVNv SB + Uby 5Y+ = 



Um = -M 2 ad(A+)ad(A_) , 

Uby = — £tad(A_)p* + ^ad(p»A_) . 

At this point, we collect the results obtained so far as follows: Let 



$ 



4>y \ I SY + 

4>b ) \ VNrSB 

then ( 3.33| ) and ( 3.37 ) can be written as a 2 x 2 matrix equation 

$ + [/$ = 

with 

Uyy Uyb 

Uby Ubb 



U 



(3.37) 

(3.38) 
(3.39) 

(3.40) 

(3.41) 

(3.42) 



The operators in this matrix are given in Eq. (3.34), ( |3.35 ), ( 3.38| ) and ( 3. 39] ). 
The perturbation equations ( 3.30| ) and (3.41) do not include the Gauss con- 
straint (p. 24), whose linearization is easily found to be 



0, 



I p*- <f>B + ad(A_) (f> Y J 



= 0. 



(3.43) 



The role of this constraint will be discussed below. 

In concluding this section, we emphasize once more, that the perturbation 
equations hold also for black holes, if these are assumed to be of essentially 
magnetic type (see Eq. (3.1)). We also would like to note that a comprehensive 
discussion of the pulsation equations for the SU(2) YM-Higgs sphaleron can be 
found in Ref. lEl. 
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4 Transformation to a hyperbolic system 

A look at the second order differential operator U shows that it is not elliptic 



and, thus, the system (3.41) of partial differ ential equations is not hyperbolic. 



With the help of the Gauss constraint ( 3.43 ) it is, however, possible to derive 



a hyperbolic system for the subspace of physical perturbations orthogonal to a 
space of pure gauge modes. This reformulation of the perturbation equations 
will turn out to be very useful for several purposes. 

We need first some notation. It is natural to introduce the following scalar 
product for LGc-valued functions on (r ,oo) C 1Z + : 

/•OO 

(ct>\^)= / (<^> dr. (4.1) 

Jra 

with the weighted measure 

dr 

For a black hole, the lower limit ro is the radius of the horizon and for a regular 
solution it is zero. The operators Uxx and U are symmetric with respect to 
this scalar product on a dense domain of L 2 -functions. This can be seen easily, 
using 

(^\p^) = (pA\^) (4.2) 

for smooth functions, which vanish at ro, and 

(0|ad(Z)^)=-(ad(c(Z))0|^) (4.3) 



for arbitrary LGc-valued functions <fi, ip, Z in L (see (2.5)) 



A "gauge mode" 4>g is by definition a perturbation of the form 

$C = - iGx , (4.4) 

where Q is the linear operator 




QX = 1 (4-5) 

V ~^ x ) 

and x is a (E)c-valued function. Note, that such variations arise if ([2.6|) is 



subjected to (T- valued) gauge transformations g — exp(— e\)- Eqn. ( |2.7D and 
( |2.§| ) show that this induces the infinitesimal transformation 

A+ ->■ A + - ad(A+)x , y/NrB-^y/NrB-i-p^x- (4.6) 

It is not surprising that the following identity holds 

UQ = Q, (4.7) 
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whence 

U$ G = . 

"Physical perturbations" $p satisfy by definition 

G$p = , 
where Q is the linear operator 

C?$=p*-(/> e + ad(A_)0y . 



(4.8) 



(4.9) 



(4.10) 



The component 4>y is assumed to have values in the subspace (2.10) of LGc 
and (fig has to be ( S )c-valued. Hence, physical perturbations are by definition 
those, for which the curly bracket in (3.43) vanishes. 

Roughly speaking, a physical perturbation is orthogonal to all gauge modes. 
More precisely, modulo boundary terms we have 



i($ P I $ 



G , 



:*p\gx) 



[G®p\x)=0, 



(4.11) 



which follows easily with Eq. (4J2) and 
The identity 

GU = 0, 



(4.12) 
which can be verified by direct calculation, is related to the Gauss constraint 

dtG$ = Q (4.13) 



in the following way: Assume Eq. (4.13) is satisfied for t = to, then the dy- 
namical equation (3.41) implies that (4.13) is satisfied for all times. Indeed, we 
conclude with (4.12) that 



dt 2 (G$) = G{d t 2 ^) = -G{U<5>) = 



(4.14) 



As a corollary we have: A solution of (3.41), which lies initially in the physical 
subspace (4.9) and satisfies initially the Gauss constraint (4.13), will satisfy the 
"strong" Gauss constraint (fOJ) for all times. For physical perturbations we can 
thus use this strong form to bring Eq. (3.41) to a hyperbolic form. After some 
manipulations, one finds 



where 



U = 



V 



{p, 2 + v}-G^ 2 G 

Vyy Vyb 
Vby Vbb 



(4.15) 



(4.16) 
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is the following (matrix-valued) potential 

V YY = H 2 K 2 + ^ 2 ad(iF||), (4.17) 

Vyb = 2{p*fi)K + + 2n,&d{p*A + ), (4.18) 

V BY = -2{p*n)K_ +2/iad(p,A_) , (4.19) 

Vbb = ^ 2 K 2 -^^, (4.20) 



with 



/' 



K 2 = -ad(A+)ad(A_), (4.21) 

K± = ±ad(A ± ) . (4.22) 



Modulo the strong Gauss constraint Q<& = 0, Eq. (3.41) is thus equivalent 
to 

<9 t 2 $ = - jp* 2 + ^}$ . (4.23) 

This system is clearly hyperbolic. We emphasize that this new system implies 
the strong Gauss constraint for all times, if it is satisfied ini tial ly: Q $ | tn = 



Q dt$ \tn = 0. The argument runs as follows: As a result of (4.8), ( 4.15 ) and 



(122), Q<& satisfies the hyperbolic equation 

d t \G<S>) = -QG^rn) = -{p*~^P* +K 2 }n 2 (g$) . (4.24) 



Uniqueness of the Cauchy problem for the hyperbolic system (4.24), with ap- 
propriate boundary conditions at ro, then implies our claim. 

We specialize now to harmonic perturbations proportional to e~ tuJt and ob- 
tain for the amplitude of <&, denoted by the same letter, the two eigenvalue 
problems 

C/$ = cj 2 $ (4.25) 

and 

|p* 2 + V j$ = w 2 $ . (4.26) 

The second equation has the form of a (vector-valued) Schrodinger equation. 

In the next section, we prove that the spectrum of U has a nonempty negative 
part (which is presumably discrete) , by constructing a smooth trial function (5$ 
for which ( <5$ | U | <5<I> ) is strictly negative. This implies that the operator p 2j rV 
has also a negative part in the spectrum. This can be seen as follows: 

If we can show that there exists a smooth function \, such that 

igs$ = gg x = {p*-^p* + K 2 }x, (4-27) 
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then we have a decomposition 

S<P = 8<S> P - iQx (4.28) 

into smooth physical and gauge components. Using also ( [1.8] ), we have 

( <5$ | U | 5$ ) = {5®p | U | £$p ) = ( 5<P P | p, 2 + V| 5$ P ) < , (4.29) 

which would imply our claim. 

Since QQ is a positive operator and since iQ6<& is smooth, we expect on 



the basis of elliptic existence and regularity theorems that ( 4.27 ) has indeed 
a smooth solution. This is one of several mathematical points which will be 
discussed in the appendix. Another issue will be, whether the operator pj 2 + V 
is essentially self-adjoint on a dense domain of smooth functions, which satisfy 
the boundary conditions implied by the physics of the problem. This will be 
analyzed in section 6 and in the appendix. 

The relation between the operators U and Q :— p* 2 + V, given explicitly in 
( 1.15J ), can be summarized (on a formal level) as follows: As a result of (1.7) 



and (4.15), both operators split relative to the decomposition of the L 2 space of 
perturbations into physical and gauge degrees of freedom, L 2 = Tip ® TCg, an <4 
their restrictions satisfy 

Q\h p = U\h p i 

U\h g - 0, 

Q\n G = Q^QUc > ■ 
The last inequality follow from 

{<S>g\Gh 2 G\<5>g) = {G<S>g\" 2 \G®g) 

for <&g G 7~Lg- In particular, the negative part of the spectra of U is contained 
in that of Q and the discrete spectra of the two operators coincide. 

5 Instability of generic EYM solutions 

We are now ready to establish the main point of this paper: For a given so- 
lution with A + = X^aes w a e a, we construct a one-parameter family of field 
configurations A( T ) + , B^ such that (8&\U\8<&) < for the variation 

,HI>=| '] = | ^ A M+l-° | , ,.-,.!, 

SNrd T B (T) \ T = 
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The families we consider are of the form 

A (t)+ = Ad(exp(rZ))|A + cos(r)+iT + sin(r)| , (5.2) 

B(r) = -tZ', (5.3) 



where T + is a real element in the subspace (2.10), satisfying 

[T + ,T_]=-2tA3||, (5.4) 

and Z is a ( E ) -valued function of r with 

lim ad(A + )Z = iT + , suppZ' c (ro,oo) . (5-5) 

r — >ro,oc 

If E is not empty such an element T+ always exists (see appendix A of ref. 
m). A function Z with the required properties can be found if 

lim w a 7^ for all a E E . (5-6) 

r — >ro,oc 

This can be seen as follows: Let {/i Q }aes be the dual basis of 2-7rE and put 

Z = ^2 Z a h a , T + = ^ T a e a (5.7) 

and 

7 _/ -T a /w a (r ) for r < r + (1 - e) , , 

Q ~ 1 -T q / Wq (oo) for r>r + (l + e) ^ 8j 



for an e > 0. Then, both conditions in ( |5.5| ) arc satisfied. 

For a regular (uncharged) solution, condition (|5.6D is fulfilled and E is not 



empty H. Thus, a family (5.2), (5.3) always exists for solitons H. 

We note some properties of the families above. For the gauge group SU(2), 
these are closely related to families studied by other authors pa. The equilib- 
rium solution is clearly obtained for r = 0. Applying a gauge transformation 
with g — exp(rZ), we obtain 

A (r)+ -> A+ cos(t) + i T+ siruV) , B {t) -> . (5.9) 



The first variations of (|5.2|) and (5.3) are 



6(f>Y =iad(A+)Z + T+, 5<t> B = -i-p*Z , (5.10) 

and these satisfy by construction the desired boundary conditions 

lim 54> Y =0, lim 506 = 0. (5.11) 

T— »r o,oo r—^r 0,00 
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(S(f>B has even compact support in (ro,oo)). Since an equilibrium solution sat- 
isfies 

p*A+\ r „ = p*A+\ 00 = 0, (5.12) 

we also have 

lim p*84>Y = . (5.13) 

r — >ro,oo 

This choice of trial functions fulfills our goal: 5$ is normalizable and 
( <5$ | U | 5$ ) is finite and turns out to be strictly negative. 



The first of these two points is simple. Since <50e in (5.10) has compact 
support, we have to check only whether 



l<5<M 2 ^ < oo 



(5.14) 



By construction, 

^ T a (w a /w a (r ) - lj e a 



for 



< r + (1 - e) , 



by = > a£E 



J^ Taiwa/wa^oo) - lj e Q for r > r + (1 - e) 



(5.15) 



a£S 



Hence, the integrand has a finite limit for r — ► ro (even for extreme black hole 
solutions) . Since iV and <S* both approach one at infinity, the integral is finite if 
A + — A + (oo) converges to zero faster than r~ x l 2 . 

The calculation of ( 5$> | U \ 5Q ) is somewhat tedious. Considerable simplifi- 
cations occur by separating a gauge mode in 5$: 



5$ = (5$ - iQZ 



with 



<5$ = 








(5.16) 



(5.17) 



We stress, that neither (5$ no r QZ are normalizable. Nevertheless, we have 
UQZ = and thus (5.16) and (5.17) give (with a slight abuse of notation) 



'8$\U\8$) 



= (5®\U\5®)+i{gZ\U6®) 

= (6$\U\6$)+i( UQZ \8$) + ( ad{ P *A+)Z , T + 

= (6$\U\6$). 



ro 

(5.18) 



The boundary term does not contribute because of Eq. ( |5.12 ). From this, we 
obtain the intermediate result 



'6$\U\6$) = 



u, 2 (T + , ad(i^||) T + ) dr* 
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M 2 (^||, A 3 



dr* 



(5.19) 



where we have used ( 2J5 ) and the property (5.4) of T + . 
Finally, we show, that the last term has a definite sign: 



IJ?(f\\ ) A 3||) dr* = 



|p*A+| 2 



2 J u 2 |J tr ||| 2 dr* 



(5.20) 



After a partial integration, we find with the unperturbed YM equation (2.23) 



|p*A+| dr* 



-z(p»A + , A+> 



M 2 (A+, a,d{if\\)A+)dr*. (5.21) 



The boundary term vanishes because of Eq. (5.12|), and since 



2|^|| |" = (^||,*[A + ,A_]-2A 3 ||) 

= (A+, ad(i£y)A+) -2(^||, A 3 ||), 

we have established the crucial result 

(8$\U\8$) = -(p*A+|j»*A + )-2(/i^||| i ^||) 

/°° r 2 i 



(5.22) 



(5.23) 



This expression is clearly finite an d str ictly negative. 

One can show, that expression ( 5.23| ) is also equal t o th e s econ d variation of 
the Schwarzschild mass for the one-parameter family (5.2), (5.3). (This is the 
way we arrived originally at the variation (5.10)). For a systematic discussion 
of the relation between variational principles for the spectra of radial pulsations 
and second variations of the total mass, we refer to ]2q |. 

In summary, we have proven (apart from technical subtleties) that static, 
spherically symmetric, asymptotically flat solutions of the EYM equations are 
unstable. More precisely, we have established: 

Theorem 1 A generic, regular solution is unstable, if the (magnetic) YM 
charge vanishes (i.e., if linv^oo A(r) is a homomorphism from LSU(2) to LG) 
and if asymptotically A + — A + (oo) ~ r~ a with a > 1/2. 

For a black hole (with horizon at th and A + = X)qgs w a e a)i the assumptions 
are somewhat more restrictive and "trivial" solutions have to be excluded. We 
call a generic, essentially magnetic solution "trivial" if either £ is empty or 
each amplitude to a is constant. These are clearly just the Reissner-Nordstrom 
solutions. 
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Theorem 2 A generic, essentially magnetic, non-trivial black hole solution is 
unstable, if ]im r - trh>oc w a 7^ for alia € S and if asymptotically A + — A + (oo) ~ 
r~ a with a > 1/2. 

We would like to stress that we were able to draw this conclusion, assuming 
only weak asymptotic conditions for the solutions. In particular, the fall-off 
condition is mild and is certainly fulfilled for the Bartnik-McKinnon and the 
related black hole solutions, as was shown rigorously in [£7j. The same is true 
for the regular solutions, which have been found numerically by H.P. Kunzlc for 
the group SU(3) |2j. (For both types, the exponent a is equal to one.) 

6 Sphaleron-like instabilities as bound states of 
a fictitious deuteron problem 

We address now the question, whether the operator p\ + V in the eigenvalue 



problem ( |4.26| ) is essentially self-adjoint on a dense domain of smooth functions, 
which satisfy the boundary conditions implied by the physics of the problem. 
That this is indeed the case, will be shown in the present section for SU(2) 
solitons. The discussion of the general case is deferred to the appendix. 



For regular SU(2) solutions, it turns out, that the eigenvalue equation ( |4.26|) 
can be interpreted as a fictitious deuteron problem for a neutron-proton po- 
tential, consisting of a central part, a tensor force and a spin-orbit coupling. 
All parts are determined by the unperturbed soliton and can be shown to be 
bounded. The corresponding Schrodinger operator is thus essentially self-adjoint 
on the subspace of smooth functions with compact support and self-adjoint on 
the Sobolev space H 2 (TZ 3 ). These facts will be used later in an analysis of the 
instabilities, implied by the existence of bound states (see the appendix). 

In order to bring the operator pi + V to a standard Schrodinger form, we 
introduce the new radial coordinate 

*>=(! <") 

in terms of which p„ = —id/ dp. Since u 2 behaves like l/p 2 near the origin, we 
separate from the potential ( |4.16| ) the singular term 

V(p) = J — + V( P ) , (6.2) 



l>- 



where 



f-\ K '' 0) 2,A ' +<0) | (6.3) 

-2iKL(0) K 2 (0) + 2 



and the remainder V is bounded. 
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For a generic soliton, the eigenvalues of J 2 are equal to j(j+l) with j — k±l, 
whereby the integer k runs through a (strictly) positive, finite set. This set 
always contains k = 1 and is uniquely determined by A3 = A311. 

In a representation in which J 2 is diagonal, J 2 / ' p 2 thus describes the central 
barriers of a finite set of partial waves. 

We now discuss in detail the equation s for the gauge group SU(2). For this 
group, only S and D waves occur in Eq. ( |4.26| ). In the variable p and with the 
parametrization 

A + = w t + = w (n + ir%) , 

A 3 = T 3 

and 



4>Y 




T3 




P 

For the potential V we find 



where 



(6.4) 



the eigenvalue equation ( 4.26 ) takes the form 



3 \ v DS v DD 



V SS = {p / '/^ + 8(pw)' + 6p 2 w 2 }+2p, 2 (l-w 2 ), (6.7) 

V SD = V2{ M '7 M + 2(^)'}-V2/i 2 (l-u. 2 ), (6.8) 

V DS = V SD , (6.9) 

Vdd = 2{n"/n-4(nw)' + 3p 2 w 2 -9/p 2 }+n 2 (l- w 2 ) (6.10) 

and a dash denotes a derivative with respect to p. 

It is amusing and helpful to note, that this coupled eigenvalue problem has 
the same form as the Schrodinger equation for the relative motion of a two-body 
proton-neutron system with the three standard terms Vc(r) (central potential), 
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Vt{t) S12 (tensor interaction) and Vz,s(r) LS (spin-orbit interaction). For total 
angular momentum J = 1 and total spin S = 1, the possible orbital angular 
momenta are L = 1 and L = 0, 2. Because of parity conservation, the P wave 
decouples from the S and D waves. The remaining equation, describing coupled 
S and D waves, reads in suitable units 




Vc V8V T \\{ Us \=f( Us \ (6U) 

V8V T V c - 2V T - 3V LS J J \ud J \u D J ■ [ ■ > 

These equations have first been derived by Rarita and Schwinger pa]. Our 
eigenvalue problem (p^a) is clearly just a special case of (|6.11) and we can, by 



identification, express the three potentials in terms of the functions N, S, w of 
the Bartnik-McKinnon solutions. 

We present numerical results elsewhere (see also Ref. p9|] ) and emphasize 
here only, that with this interpretation the mathematical nature of our eigen- 
value problem is automatically settled, because the perturbation V is completely 
harmless. We come back to this in the appendix, where we discuss also the op- 
erator corresponding to the strong Gauss constraint. 

Acknowledgments 

Special thanks go to Dieter Maison for very useful comments. Interesting dis- 
cussions with members of our theory group, especially with George Lavrelashvili 
and Michael Volkov, are gratefully acknowledged. We would also like to thank 
Markus Heusler for discussions at an earlier stage of our work. Finally, we wish 
to thank the Swiss National Science Foundation for financial support. 



Appendix 

In the main body of the text, we deferred on several occasions some of the 
mathematical subtleties to this appendix. 

A Essential self-adjointness of the effective 
hamiltonian 



Fo r bla ck holes, the operator Q = pi + V in ( 4.15|) , with the expressions ( 1.17 ) 



( [1.20 ) for the matrix- valued potential V, is effectively a standard Schrodinger 
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operator on the whole real line (see JL7|) and is thus essentially self-adjoint on 
C°° functions with compact support. (The potential V is bounded for black 
holes.) For solitons, we can use Weyl's limit point - limit circle criterion (see 
[ |30| . Sec. X.l or pi]]) for the first two terms of the operator 

Q =- t +J i +V{p) (A1) 



(see (|6.2|), (6.3)). Since V is bounded, the Rellich-Kato theorem implies, that 



the domains of (essential) self-adjointness are not changed by this additive term. 
Another method which will be used later, is to lift Q to a Schrodinger op- 
erator Hq on TZ 3 and to use powerful results for this kind of operators. In Sec. 
6 we showed how this can be achieved, if the gauge group is SU(2): Hq can 
then be chosen to be of the standard form for a deuteron problem. This oper- 
ator is essentially self-adjoint on C (1Z 3 ) ® C 4 and self-adjoint on the Sobolev 
space K 2 (7Z 3 ) <g>C 4 (see, e.g., Q, Sec. X.2). Restricting these domains to the 
subspace of S and D waves, provides the domains we are interested in for the 
original operator Q. For instance, Q is essentially self-adjoint on 

V(Q) = { (u s ,u D ) | u s G C§°[0,oo),u s (0) = 0;u D e Cg°(0,oo) } • (A2) 

Although, we have not yet generalized this construction to arbitrary gauge 
groups, the generalization of T>(Q) is obvious: The S waves have to be restricted 
as in ( |A2| ) and the higher waves have to lie in C^°(0, oo). We also note at this 
point that the variation ( [5.1C| ) lies in the domain of definition of the self-adjoint 
extension of (Q,T)(Q)). 

If we would restrict the S waves also to C^°(0,oo), the operator Q would 
not be essentially self-adjoint. For each S wave sector, it would actually have 
a one-parameter family of self-adjoint extensions. The self-adjoint extension, 
given above, is just the Friedrichs extension, and one can show that it is the 



only one which is compatible with the strong Gauss constraint (4.9). 

The existence and smoothness problems in connection with Eq. (4.27) can 
also be solved by lifting the equation to TZ 3 and using standard existence and 
regularity theorems for elliptic operators on TZ 3 . (The details can easily be 
worked out for G = SU(2).) 

B Spectral properties and unstable perturba- 
tions 

In Sec. 4 it was shown that the perturbation equations for even parity pertur- 
bations are equivalent to the hyperbolic system 

3 t 2 $ = -Q$. (Bl) 
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We recall also that these equations imply the propagation of the strong Gauss 
constraint. As a main point of this paper we proved that the self-adjoint opera- 
tor Q, restricted to the subspace of physical states, satisfying the strong Gauss 
constraint, has a non-empty negative spectral part. This fact implies, of course, 



that there are unstable Hilbert space solutions of (Bl). We just have to choose 
the initial data $o such that Eq(— oo,0) $o 7^ 0, where Eq(-) denotes the pro- 
jection valued measure belonging to Q (see below). It is even possible to choose 
$o € C^°, because the smooth functions with compact support are dense in the 
Hilbert space L 2 . 

The question now arises, whether such a Hilbert space solution is even a 



(classical) solution of the system of partial differential equations flBl| ) , in other 
words, whether the Hilbert space solutions with $o € C^° are automatically 
smooth. For black holes, the positive answer to this question is contained in a 
paper by Wald |32]. His analysis does, however, not directly apply to solitons, 
because he assumed, that space is a complete Riemannian manifold. 

A direct attack of the problem on the half-line (0, oo) is again difficult. Once 
more, a way out is lifting the problem to ll 3 , 

d t 2 <P = -H Q $ , (B2) 

where the analysis of |32J applies. We showed earlier, how this can be done for 
SU(2). Since the required smoothness properties certainly do not depend on 
the gauge group, it is not worthwhile to elaborate further on this. We would 
like, however, to present here a simplified version of Wald's argument. 



Consider a hyperbolic system on 1Z x W l of the form (B2) with a smooth 
elliptic operator A (instead of Hq), which is essentially self-adjoint on C§°(7?.™). 
For systems of this kind, a lot is known about the Cauchy problem (a standard 
reference is p3|). In particular, one knows (see Theorem 23.2.2 in Vol. Ill of 
p3|), that the Cauchy problem with initial data <I>o, ^o G &^'{W l ) has a unique 
solution in C°°(K xK n ), which for any fixed time t is in C^(K n ). This smooth 
solution must agree with the Hilbert space solution of the Cauchy problem 
because the latter is also unique. 

Let us now assume that the spectrum o~(A) of A has a non-empty intersection 
<r{A)- with (—oo,0). We also assume that <j{A) is bounded from below (this is 
the case for A = Hq). The Hilbert space solution of the Cauchy problem can 
easily be expressed in terms of the projection valued measure E(-) belonging to 
A. It suffices, for what follows, to take as initial data $|t=o = $o,d t $\t=o = 0. 
Then, the corresponding Hilbert space solution of 

$t = -A$ t (B3) 



is 



$ t = E({0})$ + / cos(tVA) dE(A)$ 

J(O.oo) 
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I cosh(iV^A) dE(A)$ , (B4) 

Jer(A)_ 



l<r{A)_ 

as can easily be verified. Note in particular, that &t=o = E(7?.)$o = "£o, as 
required. 

With standard rules (see, e.g., [p4J Chap. 13), we obtain from this 



;$ |$ f > - ||E({0})$ || 2 + I cos(*\/A) d/i* (A) 

J(O.oo') 

cosh(i\/— A) <i/j$ (A) (B5) 



*(A). 



and 



||$ t || 2 >/ cosh(tV^A) d//*„(A), (B6) 

where /z$ is the finite measure ( $o | E(-)$q ) on 72., whose support is contained 
in <j(A). As emphasized above, we can choose $o such that (supp^$ ) n tr(A)_ 



is non-empty. Then ( |B5| ) and (B6) imply that both quantities on the left di- 
verge exponentially. This exponential grows translates to an average exponential 
grows of the classical solution of the hyperbolic system for smooth initial data 
with compact support. 

These considerations conclude our instability proof. 
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